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Abstract 


We have devised a technique that makes it possible to obtain energy estimates for 
initial-boundary value problems for nonlinear conservation laws. The two major tools to 
achieve the energy estimates are a certain splitting of the flux vector derivative f(u) X) and 
a structural hypothesis, referred to as a cone condition, on the flux vector f{u). These 
hypotheses are fulfilled for many equations that occur in practice, such as the Euler 
equations of gas dynamics. It should be noted that the energy estimates are obtained 
without any assumptions on the gradient of the solution u. The results extend to weak 
solutions that are obtained as pointwise limits of vanishing viscosity solutions. As a by- 
product we obtain explicit expressions for the entropy function and the entropy flux of 
symmetrizable systems of conservation laws. Under certain circumstances the proposed 
technique can be applied repeatedly so as to yield estimates in the maximum norm. 




1 Introduction 


Most of the existing theory for nonlinear conservation laws is concerned with the initial 
value problem. The basic tenet of this paper is to devise a technique that makes it possible 
to obtain energy estimates for the initial-boundary value problem. The key to obtaining 
an energy estimate lies in a certain splitting of the flux vector derivative /(u) x . Based 
on this splitting one can obtain energy estimates without assumptions on the gradient 
u x , provided the flux vector / satisfies certain structural hypotheses. It should be noted 
that these hypotheses are fulfilled for many equations that occur in practice, such as the 
Euler equations of gas dynamics. In certain cases the proposed technique can be applied 
repeatedly so as to yield estimates in the maximum norm. 

We now give a brief presentation of the contents of this paper. Section (2) deals 
with the scalar problem in one space dimension. First, we consider flux functions of 
the form f(u) = u\ for which it is very easy to derive an energy estimate. Next we 
assume that f(u) can be expanded in a convergent Taylor series and apply the technique 
to each individual monomial vP . This procedure leads to a certain differential equation. 
Conversely, the solution to this differential equation will achieve the correct splitting of 
/(u) x , even if / is non-analytic. To obtain an energy estimate, the “dangerous” boundary 
terms must be eliminated by means of the characteristic boundary conditions. This can 
not be done in general, however. We therefore propose a so called cone condition to single 
out the relevant fluxes. In section (3) we repeatedly apply the ideas developed in the 
previous section. The main result is a maximum norm estimate for the scalar initial- 
boundary value problem. The principles of section (2) can be generalized to systems in a 
straightforward manner. This is the major issue of section (4). To begin with, we confine 
ourselves to symmetric hyperbolic systems, i. e., systems where the Jacobian satisfies 
fiT _ j-i The cone condition is generalized to systems, and we are thus able to derive an 
energy estimate. Next, we require that the system be symmetrizable by means of a change 
of variables u = u(v). The Euler equations belong to this latter category. Applying the 
splitting to the time derivative as well, we are able conclude that v satisfies a generalized 
energy estimate. We use the term generalized energy estimate, because the norm will, 
in general, depend on u(v) (but not on the gradient of u). As a by-product we get an 
explicit expression of an entropy function t/(u) and the corresponding entropy flux f(u). 
Section (5) is concerned with generalizations of the preceding sections to several space 
dimensions. Finally, in section (6) we prove energy estimates for weak solutions that can 
be obtained as limits of vanishing viscosity solutions. 
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2 Scalar Conservation Laws in One Space Dimen- 
sion 


Consider the conservation law 


u t + fx — o, X 6 (0, 1) t> 0 
u(x, 0) = ip{x) . 


(1) 


At the boundaries x = 0, 1 we prescribe data ip(t) for the ingoing characteristics, which 
are determined by the sign of /'(u(i, t)), i = 0, 1. The flux f(u ) 6 R is assumed to be 
a continuously differentiable function of u € R. It is well known that eq. (1) allows an 
energy estimate when f(u) = u 2 . One way to see this is to make use of the splitting 


(” ! ), = \ 1(” 2 ), + uu * 

which is valid for smooth functions u. If we use this splitting in eq. (1) and apply the 
energy method we obtain 


^IM| 2 = -^ [ u ( u2 ) r + « 2 «*] dx = ~\j Q (u 3 )^dx = ^ (u(0,f) 3 - u(l,f) 3 ) , 

where || ■ || denotes the usual norm associated with the scalar product 


(u, v) = I uvdx . 

Jo 

If u(0,<) > 0, then we prescribe u(0, t) = Similarly, if u(l,/) < 0, then we set 

u(l,<) = We then obtain an energy estimate. 

More generally, any monomial f(u) = u J satisfies the identity 




(2) 


Suppose for the moment that f(u) is analytic. Then there exists a convergent Taylor 
series such that 

OO 

/(u) = £ CjU 3 . 

j=0 

Since the convergence is uniform we get 


using eq. (2). Some simple calculations show that 

,22. j 1 ru 

f f(v)dv. 

j=o J + 1 u 


( 3 ) 


(4) 
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Define 


( 5 ) 


F(u) = - r f(v)dv= r f(9u)d6. 

u Jo Jo 


Eqs. (3), (4), (5) imply that 

f, = U-F)r + -(f-F)«,. (6) 

u 

From definition (5) it follows immediately that F satisfies the following differential equa- 
tion 

F'u = -F + f m 

F(0) = /(0). 

It will henceforth be assumed that /( 0) = 0. This is no restriction, since we can always 
replace f(u) by f(u) — /( 0) without affecting eq. (1). Furthermore, most examples that 
occur in practice satisfy /( 0) = 0. Combining eqs. (6) and (7) gives 

f x = (F'u) x + F , u x , (8) 

which we will refer to as the canonical splitting. 

Conversely, let F(u) be the solution to eq. (7), where it is no longer required that f(u ) 
be an analytic function. It is clear that F(u) is given by eq. (5). Hence, 

f: = (f-F) x + F x = (F'u) x + F'u x , 

where the second equality follows from eq. (7) and the chain rule. We have thus established 
the canonical splitting for all / such that eq. (7) has a C'-solution. We shall henceforth 
assume that / is differentiable. We point out that if / is a linear function of u, i. e., 
f(u) = au, then F(u) = (1/2 )au and thus F'(u ) = (1/2 )a. Consequently, 

, 1, . 1 
Jr — 2 T 2 i 

which is the usual skew-symmetric form of a linear hyperbolic PDE. 

The energy method applied to eq. (1) yields 

4n| 2 = -2(«, fr) = — 2(u, (F'u),) - 2 (u, F'u x ) = -2uF'u\l , 
at 

where the second equality follows from the canonical splitting (8). To obtain the actual 
estimate we must analyze 


iF'u = uf — ! f(v)dv= f f'(v)vdv. 
Jo Jo 


From the last expression it is clear that, in general, knowledge about the sign of /'(u), 
upon which the characteristic boundary conditions are based, has no bearing on the sign 
of uF'u. We must thus limit ourselves to flux functions / such that the sign of f'(u) will 
determine the sign of 

f f(v)vdv = uF'u , (9) 

Jo 

if we are to obtain an energy estimate. 
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Definition 2.1 A function f : Q — ► R, fi C R, is said to satisfy a cone condition if 

sgn(f'(u)) = sgn i^J f'{v)vd\^j uffl. 


This cone condition is satisfied if 


sgn(/'(u)) = sgn(u), u6(l. 


(10) 


Suppose that f'(u ) > 0 at the lower boundary x = 0. We then have an ingoing character- 
istic, which implies a boundary condition at x = 0. Now, eq. (10) implies u > 0, whence 
0 < v < u. But then (10) implies that f'(v) > 0, that is, f'(v)v > 0. Consequently, 


uF'u = f f'(v)vdv > 0, 
Jo 


which shows that the boundary energy has the needed sign. Similarly, if f'(u) < 0, i. e., 
there is an outgoing characteristic, then u < 0 from (10), which again leads to f'(v)v > 0. 
But the upper limit of the integral in eq. (9) is negative, and so the entire expression is 
negative. Below we show the graph of f'(u) for a function / satisfying the condition (10). 



It should be noted that the condition (10) holds for all monomials of even order; in 
particular, it is true for Burgers’ equation. Furthermore, it is also true for the non-convex 
Buckley-Leverett equation for u € [0, 1], which is a natural restriction since u describes 
the saturation of water in a two phase flow. 

Another sufficient condition for the cone condition to be true is furnished by 

sgn(/'(u)) > 0 (or sgn(/'(u)) < 0) u 6 Cl . (11) 

For, if f(u) > 0, then the integrand of eq. (9) satisfies sgn(/'(u)t;) = sgn(v). Hence, 



sgn(/'(u)) = sgn Qf f'(v)vdvj , 

which shows that / satisfies a cone condition. Monomials of u of odd order, linear fluxes 
in particular, belong to this second category. 

Before proceeding we note that the boundary data r/>(f) may be restricted; for example, 
the cone condition (2.1) requires that /' and uF'u have the same sign. This is the case 
for Burger’s equation (/'(it) = u,uF'u = (l/3)u 3 ) in which case the condition (10) holds. 
If we have an ingoing characteristic at x = 0, then clearly u > 0 and we may thus only 
prescribe data ip(t) > 0. From now on it will always be assumed that the boundary data 
is compatible with the cone condition. 

Proposition 2.1 Let u be a smooth solution of the initial-boundary value problem (1), 
where it is assumed that f is differentiable and satisfies a cone condition (2.1). Then u 
satisfies an energy estimate. 


Proof: 

By means of the canonical splitting we get 


f t \\uff = -2uF'u\ 


1 

0 5 


where 


uF'u = / f'(v)vdv . 
Jo 


Since / satisfies a cone condition it follows that the sign of uF'u is uniquely determined 
by that of /'(u). For convenience, we assume that there is an ingoing characteristic at 
x = 0 and an outgoing one at x = 1, i. e., /'(u(0, t )), /'(u(l, t)) > 0. Hence, 


d f&tt) 

5 l|u|| 2 <2jf /'(»)«*>, 


where the data tp(t) at x = 0 is such that the integral in the right member is nonnegative 
(because of compatibility with the cone condition). Integration with respect to time yields 

n t P(t) 

f (y)vdvdr , 

J 

which proves the proposition. D 

Remark: If there are no ingoing characteristics we get an energy estimate in terms of the 
initial data alone. Also, if there were an ingoing characteristic at x — 1 (/(u(l,t)) < 0), 
then there would be a second boundary term in the energy estimate. In the case of a 
linear equation (/ = au) we get uF'u = (l/2)ait 2 , and we thus recover the familiar form 
of the energy estimate. It should be emphasized that no gradients appear in the estimate; 
if f'(v) oc v p-2 , then there is an energy estimate for all time if xp(t) G L p ( 0, oo). □ 
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3 Sequences of Conservation Laws 


The canonical splitting (8) implies that the original conservation law (1) can be formulated 
as 

u t + ( F'u) x + F'u x = 0, 

where F satisfies eq. (7). Multiplication by u and the product rule yield 
where 

fi — uF'u = f f'(v)vdv. 

Jo 

This is a new conservation law for (1/2 )u 2 associated with the flux /i = uF'u. But then 
we can split this new flux in canonical form 


where F\ satisfies 

F{u = -F l +f 1 

Fi(0) = 0 . 

Hence, multiplication of ((l/2)u 2 ) ( + (/i) r = 0 by u implies 



(3« 3 ) t + (/2)*-°, 

where 

ftt 


/ 2 = uF[u = / f[(v)vdv. 
Jo 

But 

fl = ( uF'u y = ( uf - uF)' = f - F - F'u + f'u = 

whence 

h = / f'{v)v 2 dv . 
Jo 


The process can be repeated. A simple induction argument shows that 

+ (/,)* = 0, p = 0,l,... , (12) 

where 

f P = [ f(v)v p dv, p = 0 , 1 ,.... 

Jo 

We have thus established a sequence of conservation laws. In the following we shall restrict 
ourselves to the subsequence where p is odd, i. e., 

+(/ 2 P -i) I = 0, f 2p -i = f'(v)v 2r> -'dv. (13) 
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The conditions (10) and (11) guarantee that 

sgn(/ 2p _i) = sgn(/') , p= 1,2,.... 

Consequently, integration of eq. (13) over (0, 1) x (0,<) yields 

ft f 4 >(r) 

< llvll S + 2 p/ 0 / 0 f'(v)v ip ~ l dv<tr , p = 1,2,..., 

where we have assumed /'(u(0,i)), / ; (u(l, t)) > 0 for convenience. The elementary in- 
equality i i 

(x + y)n < X* + y» , x, y > 0 , 


( 14 ) 


implies 


( ft f^(r) \ 2 P 

|«(-, 0 l| 2 p < IM Up + ( 2 P J 0 J 0 f'(v)v 2p -'dvdTj , 


(15) 


where the last term of the right hand member is nonnegative since we have required that 
the boundary data be compatible with the cone condition. But 


and 


( 2p X Jo If'iv'jlv^-'dvdT 


2 P 


< 


max 

0<t'<|V'( T )l 

V °< r <* 


l/»l £ mr)? r ir 


\ 


2 P 


Define the maximum norm of the boundary data p G L°°(0, oo) as follows 


ilV’IU = SUp |0(T)|. 

0 <T <00 


(16) 


Then 


/ 


Wf> 


( v)v 2p x dvdr j < 1 1 0 1 1 


\ 


max |/'(t>)|* 

V o<r<< ) 


2p 


Let t be fixed and let p — ► oo. We get 

rt t M>( r) 


limsup(2p / / f'{v)v 2p l dvdr J < ||V’|I°°- 

p— >00 \ JO Jo ) 
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Similarly, it can be shown that 


lim sup ||(^|| 2 p = IHloo 

p— ► oo 

for initial data € L°°(0, 1). From eq. (15) it then follows that 

IM-,OI!oo < IMloo + IHloo. 

This estimate is valid for all t > 0, and we have thus arrived at a maximum principle for 
conservation laws. The result is summarized in the following 

Proposition 3.1 Suppose that there is a smooth solution u(x,t) to the initial-boundary 
value problem (1) fort > 0 and that f satisfies either of the conditions (10), (11). Then 
u(x,t) satisfies the maximum norm estimate 

IW-.OIIco < IMloo + IHloo. 

Remark: Lax has shown in [5] that any weak solution with compact support of the 
initial value problem (1) satisfies ||u(-,< 1 )|| 00 < ||u(-,* 0 )||oo for t x > t 0 , i. e., ||u(-,/)|| 00 is 
a nonincreasing function of time. □ 


4 Systems of Conservation Laws 


We now proceed to the case where it 6 R d is a vector of unknowns and where / = /(it) € 
R d is a vector valued function of it; each component it, is a function of x,t 6 R. We 
consider 

u t + fx = 0, i 6 (0, 1) t > 0 

u(x, 0) = (p(x) (17) 

uj(i,t) = i = 0, 1 , 

where uji and uo are the ingoing and outgoing characteristic variables. To begin with, it 
will be assumed that eq. (1) is symmetric hyperbolic, i. e., f' T = /', where 


/'(«) 


(flu, ... fxuA 
\ fd u l • • * fduj / 


is the Jacobian matrix of /. Thus, there exists an orthogonal matrix Q(u ) such that 


Q T (u)f(u)Q(u) = A(it) 


is diagonal. The characteristic variables are defined as 

u> = Q t {u)u . 


(18) 
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Let 


Q(u) = ( Qi{u ) Qo(u)) , 


where Qi contains the eigenvectors corresponding to ingoing characteristics; Qo contains 
the remaining eigenvectors, i. e., those corresponding to outgoing characteristics and those 
whose eigenvalue is zero. We now define u>i = Qju and u>o = Qo u ■ I n general, the rank 
di of Qi € R rfxrfl and the rank d 2 of Qo G R ix<i2 , d\ -f d 2 = d, will depend upon the 
solution u. 


The key to obtaining an energy estimate for the scalar conservation law lies in the 
canonical splitting (8), which in turn relies upon eq. (7). For systems of conservation laws 
we thus consider 


F' u — ~F + / piT _ pi 

m = m 


(19) 


where u, /(u), F(u) are vectors; F' is the Jacobian matrix of F. As in the scalar case, we 
require that /( 0) = 0. Equation (19) is the inhomogeneous Euler’s differential equation 
(not to be confused with the Euler equations of gas dynamics). Suppose that there is 
a solution to (19). The symmetry condition F' T = F' implies that Euler’s differential 
equation can be written as / = (u T F)\ i. e., 


/' = (u t F)" , 


which in turn shows that f lT = f. Hence, it is necessary that the Jacobian matrix f be 
symmetric. Conversely, suppose that f' T = /', then it is easy to verify that 

F{u)= P f{du)d0 (20) 

Jo 

solves Euler’s differential equation (19). Furthermore, f' T = f clearly implies F' T = F'. 
Thus, eq. (19), subject to the constraint F lT = F', has a solution (20) iff f lT = f. For 
details on how to solve (19) we refer to [4]. Note the complete similarity between eqs. (5) 
and (20). With F given by (20) we achieve the canonical splitting 


f x = {f-F) x + F x = {F'u) x + F'u x . 


Consequently, the solution u of eq. (17) satisfies 

i||u|| 2 = -2 = — 2(u, (F'u) r ) — 2(u, F'u x ) . 

Partial integration and F' r = F' imply that 

^IMI 2 = —2u t F'u\q , (21) 

which, formally, is an energy estimate. To get a true energy estimate we must make sure 
that the boundary terms 

u t F'u= f u T f'{0u)dud0 
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have the correct sign by applying the characteristic boundary conditions. As was seen 
in the scalar case, this will not be true for arbitrary flux functions. This remark also 
pertains to the vector valued case. Some kind of cone condition is needed. We generalize 
definition (2.1) to 

Definition 4.1 A vector valued function f : fi — » C R 1 is said to satisfy a cone 

condition if 

[ u T f'{9u)8ud6 > (<) 0 u£M, 

Jo 

where M Q SI is a submanifold such that 

u T fu > (<)0. 

We note that the cone condition is fulfilled if there are constants Cq > 0 and C\ > 0 
such that 

cif < F' < cof . (22) 

Assuming that this condition is true, then, by eq. (21) and the definition (18), we obtain 

— ||u|| 2 < 2 cou? t Au>(0, t) — 2cic*; r Au;(l, t) . 

But 

w^'Ao>(0, t ) — u>+A + u> + (0, t ) + loq Aou>q(0, t ) T u>^A_u>_ (0, t ) , 

where A+, Ao, A_ contain the positive, zero, and negative eigenvalues of f for a given value 
of u. At x = 0 the positive eigenvalues correspond to ingoing characteristics; the remaining 
ones are treated as outgoing characteristics. With A / = A + , A o = diag(A 0 — A_), 
w/ = u? + , and u>o = (u; 0 u>~) it follows that 

w t Aoj(0 ,t) = wjA/w/(0,<) - cjjA o ujo(0,f) . 

Similarly, at x = 1 we have 

— w t Au>(1, t) — loJ A/w/(l, t) — u>oAowo(l,f) , 

where A / = — A_, Ao = diag( A 0 A+), uj = w_, and a>o = (a?o w+)- It should be noted 
that A/ > 0 whereas Ao > 0. The characteristic boundary conditions (17) yield 

u>J A }u>i — wqAowo £; 2 u>q AjS — A o) ojo T 2'0^'A/V’ • 

The first term of the right hand member is non-positive if |5(u)| is small enough. In order 
to control the second term a bound on A /(u) is needed. The characteristic boundary 
conditions can be expressed in terms of u as 

u = Pu + Q!rP, (23) 
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( 24 ) 


where the projection P(u) is given by 

P(u) = Qi(u)S(u)Qo{u) T + Qo{u)Qo{u) T - 


Hence, 


A/(u) = A i{Pu + Qixp) . 

We, therefore, make the following 

Assumption 4.1 There exists a function A; : fi — + R d '* d ' , fi C R dl , such that 

A i(Pu + Qrt(t)) < K[(g(t)) 

holds at x — 0 and x = 1 for some vector g(t ) € R d ' , which is independent of u. 


The results can now be summarized in 

Proposition 4.1 Let. u be a smooth solution of the initial-boundary value problem (17), 
where it is assumed that f is differentiable and that condition (22) holds. If assumption 
(4.1) is true, then u satisfies an energy estimate. 


Remark: 

(i) Assumption (4.1) states that the signal speed of the ingoing characteristics can be 
estimated at the boundary independently of u. 

(ii) Assumption (4.1) is superfluous in the case of homogeneous boundary data (t/> = 0), 
since we obviously have xp T Airf = 0. 

(iii) If all of the characteristics are ingoing, then A / = A , Qi = Q, Qo = 0, which implies 
p — 0 Hence 

Ai(Pu + Qnl>(t)) = A(Q<K0)- 

Furthermore, since Q is orthogonal, there exists a vector g(t ) such that ip(t) 
Q T g(t ), i. e., Q4>(t) = g(t), which shows that 

A I (Pu + Q,xl>(t)) = A(g(t)). 

Thus, assumption (4.1) is true. D 

We point out that, in general, the boundary data must be given in characteristic form, 
even if assumption (4.1) holds, if we are to obtain an energy estimate. This is true since 

xf T Ai(u)xf < 4 > t A i{g)4> < p{Ai(g))f) T ip , 
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where ij> = Qjg for some data vector g € since Q/ has full rank. Consequently, 

ip T Ai(u)ti> < p(hi{g))g T QiQjg ■ 

Clearly, the right hand member will depend on u through QiQj. An exception to this 
is when Qj = Q, in which case we have QiQj = QQ T = I because of orthogonality. 
If, however, Q is independent of u, then of course we get an energy estimate in terms 
of the non-characteristic data g. In this context it should be noted that the coefficient 
matrices of the primitive Euler equations can be symmetrized such that the resulting Q 
is independent of u [1, 8]. 

Condition (22) can be strengthened in certain cases. Suppose that the flux vector / 
satisfies Euler’s differential equation 


f'u = pf , 


the solution of which can be written as [4] 


(25) 


/(u) = ufr(^,...,^) . 

\Ui Ui/ 


We then have 


F = f f(9u)d0 = - I' f'(eu)0udO = - I" %-BdO . 
JO P Jo p Jo do 

Integration by parts yields 

1 . W 1 1 . 1 


F =-f f f(0u)d0 = -f--F , 

P P Jo P P 


i. e., 


which strengthens condition (22) to 


F = — 
p+ V 


F' = /' . 

p+ r 


(26) 


The flux vector of the conservative Euler equations of gas dynamics satisfies eq. (25) with 
P — 1. Unfortunately, the Euler equations are not symmetric, i. e., f rT ^ /'. 

Example: We consider eq. (17) where 


Hence, 




/'(«) 
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Obviously, f' T = /' and f'u = 2/. Thus, the cone condition is fulfilled. The energy 
method yields 

^IMI ! = -§« 7 7'"I5. < 27 ) 

where we have used the canonical splitting. The eigenvalues and eigenvectors of f'(u) are 
given by 

Ai(u) = tt! - U 2 n _ ( \ _ J_ / 1 M 

A 2 (u) = u \ + u ? ’ Q ~\ q 1 92 ) _ ^V-1 I/’ 

Here we have an example where Q is independent of u. Consequently, the characteristic 
variables are 

n T 1 /u, — u 2 \ 

= Q u = — 7 = 

\/2 \ U 1 + u 2 / 

Next, we verify that assumption (4.1) holds. For brevity we confine ourselves to the lower 
boundary x = 0; the upper boundary x = 1 is treated in a similar fashion. Assume that 
Ai(u) < 0 and A 2 (u) > 0. Hence, 

A/(u) = A 2 (u) = V2lo 2 = V2Sui + y/2tl> = S\ x + V2rp < V2xf> 


if S > 0 . This implies an energy estimate provided 5 > 0 is small enough. □ 

Remark: It is not necessary that S > 0 in the previous example since 

(ui — U2) 3 + + U2) ] ^ 

(u, - u 2 ) 3 + (5 0 (tt, - ti 2 ) + v/2t^) 3 ] < V2i!> 3 

for | S | small enough. D 

We now consider eq. (17) when we no longer require that f' T = /'. Before proceeding 
we recall the definition of an entropy function. 



Definition 4.2 A convex scalar function U{u ) is said to be an entropy function for 
eq. (17) ifU(u) satisfies 

U ,T f = F ,t (28) 

for some scalar function F)u) called the entropy flux. The prime denotes differentiation 
with respect to u. 

Introducing a new dependent variable v and setting u = u(u), f(v) = f(u(v)), we get 

u t + f x = u ( + f x = u'v t + f'v x = 0 . (29) 

Equation (17) is called symmetrizable if u' = u v and /' = f v of eq. (29) are symmetric 
matrices, and if v! is positive definite. In [3, 2] it is shown that eq. (17) is symmetrizable 
iff there exists an entropy function U(u). Using the new variable v we can thus apply 
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the technique used for the symmetric hyperbolic case. The only difference is that v t is 
preceded by a symmetric positive definite coefficient matrix u', which implies that the 
change of variables u = u(u) is well defined. It is therefore natural to apply the canonical 
splitting to the time derivative as well. Thus, let 0(v), F(v ) € R d be the solutions to 

O'v = -0 + u 0 ,T = O' 

F'v = -F + f F' t = F' ( 30 ) 

where O' = 0 V , F' = F v . Hence, 

U(v) = j\{0v)d0 F(v) = jf 1 f(0v)dO . (31) 

Applying the canonical splitting to u t + f x in time and space yields 

(U'v) t + U'v t + ( F'v) x + F'v x = 0 , (32) 

which after scalar multiplication by v implies 

(y T U'v) t + (y T F'v) x = 0. (33) 

Integration of eq. (33) with respect to x yields 


-Mb = ~v T Fv\l 


where 

IMIfr — / v T U'vdx 

J 0 


defines a norm since it follows from eq. (31) that O' is positive definite. In general, O' will 
depend on v. Assume that f(v) satisfies f'v = pf for somep > 0. Then F' = (l/(p+l))/', 
whence 



P + 


j -v T fv\l • 


Let A denote the eigenvalues of /'. Arguing exactly as in the symmetric hyperbolic 
one obtains 


d 2 

J t \\ v \\l ^ — jTj- (V>o A/V»o(0,0 + V’fA/V’i (!,<)) . 


case 

(34) 


This is not an energy estimate in the usual sense because U f and A / depend on v. However, 
since are both positive definite we shall call inequality (34) a generalized energy 

estimate. We have thus shown 


Proposition 4.2 Suppose that eq. (17) can be symmetrized by a change of variables u = 
u(v) such that f(v) = /(u(v)) satisfies f?v = pf . Then any smooth solution v satisfies a 
generalized energy estimate (34 ) . 
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Remark: Harten [3] has shown that the Euler equations of gas dynamics can be sym- 
metrized by a change of variables u = u(v) such that f'v = pf. Thus, according to 
proposition (4.2) a generalized energy estimate is obtained. □ 

Equation (33) shows that eq. (17) provides an additional conservation law. It is 
therefore tempting to regard v T U'v and v T F'v as entropy and entropy flux. Indeed, we 
will show that v T U'v and v T F'v fulfill the conditions of definition (4.2). To this end we 
define the scalar functions 

U(u) = v T U'v = v T u — v T U 

F{u) = v T F'v = v T f -v T F ’ 1 ’ 

where the last equality is a consequence of eq. (30). The Jacobians U'{v ) and F'(v) are 
symmetric by construction. Furthermore, v' = v u is symmetric positive definite since it is 
the inverse of the symmetric positive definite matrix u' = u v . These symmetry relations 
and eq. (30) imply that 

U'{u) = v , F' t (u ) = v T f'(u ) . 

Also, U"(u) = v'(u) > 0, i. e., U(u) is a convex function satisfying 

U’ T f = f'T ? 


whence U(u) defined by eq. (35) is an entropy function associated with the entropy flux 
F(u). We have thus arrived at an explicit expression for the entropy /entropy flux. We 
summarize the results in 


Proposition 4.3 Let eq. (17) be a symmetrizable system of conservation laws. Then an 
entropy function U (u) and the corresponding entropy flux F(u) are given by 

U(u) = I' v T u\0v)0vd0 E(u)= I' v T f'(0v)6vde , (36) 

Jo Jo 

where v = v(u) is the change of variables that symmetrizes eq. (17). 


Remark: In the symmetric hyperbolic case we have u = v, whence U = (\/2)v. This in 
turn implies U(u) = (l/2)|u| 2 in agreement with our previous results. □ 

In the scalar case we obtained a sequence of conservation laws (13). Define the scalar 
functions 

U r (u) = ^~u 2p , p= 1, 

2p 

Hence, U' p = u 2p_1 . From eq. (13) it follows that f' 2p _ 1 = /'u 2p_1 and thus 

u' P r = n P . i, p=i,... . 
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Consequently, {U p , f 2p -i} is a an entropy /entropy flux sequence. This sequence allowed 
us to derive an estimate in the maximum norm. Assume that there exists a corresponding 
sequence in the vector valued case (17), i. e., 


(% + (% = <>, w v i r = F% P = i,..., 

where the entropies satisfy c p |u| 2p < U p (u) < C' p |u| 2p for some constants c p , C p such that 

limsup (— ) P = 1 . 

P-* 0 o \ Cp J 

We note that these hypotheses hold in the scalar case. Also, in the symmetric hyperbolic 
case we have U\(u) — U(u) = (l/2)|w| 2 . Integration over (— 00 , 00 ) x (0,/) yields 


/ oo roo 

U p (u(x,t))dx — I U p (tfi(x))dx , 

-00 7—00 


where we have assumed that u has compact support. The hypotheses on U p imply 

IM-,0lb,< gV'lMI,. 

Finally, taking the limsup of both sides of this inequality leads to the maximum norm 
estimate 

IM-.OIloo < IMIoo- 

Thus, the existence of an entropy sequence implies a maximum norm estimate for the 
initial value problem (17), provided the entropies satisfy c p |u| 2p < U P (u) < C p |u| 2p . 


Remark: Lax [5, 6] has constructed a sequence of entropy functions satisfying certain 
growth conditions for 2 x 2-systems. This sequence is then used to prove that the maximum 
norm of the Riemann invariants is a nonincreaing function of time. □ 


5 Conservation Laws in Several Space Dimensions 


We begin by considering scalar conservation laws 

n 

+ H (/»)*, =° x € C R n OO 
» 

u(x, 0) = y(x) . 

At the boundary T = dfl we prescribe the ingoing characteristics using data i/>(x, t). It will 
be assumed that the boundary T is almost regular as defined in [7]. Thus, the divergence 
theorem can be applied to fi. Let v(x) denote the outward unit normal at x £ T. The 
characteristic is said to be ingoing at x € T if 

n 

v T (x)f(u(x,t)) = ^^(x)/,'(u(x,f)) < 0. (38) 

t = l 

The cone condition (2.1) is replaced by 
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Definition 5.1 A function f : 0 U — ► R, fl u C R, satisfies an n-dimensional cone condi- 
tion if 

sgn (v T f\u)) — sgn v T f'(v)vdvj 
for all unit vectors v G R n . 


This cone condition is fulfilled if we require 

sgn(i /T f'{u)) = sgn (u T f'(v)) , 0<n<uoru<v<0. (39) 

Remark: It should be noted that the generalization of condition (10) to 

sgn(//(u)) = sgn(u), i = 1, . . . , n , u 6 (40) 

would, in general, be insufficient to ensure the cone condition (5.1), except for the special 
case of n = 1 when condition (40) implies (39). Similarly, the generalization of (11) to 

sgn(/, •(«))>(<) 0, i = 1, . . . ,n , u e flu (41) 

will, in general, fail to imply cone condition (5.1). In the special case where 

fi(u) = c x g(u ) , i = l,...,n 

condition (39) follows if g{u) satisfies either condition (10) or (11). O 

Partition the boundary asr = F/Uro, where 

Vi = {x £ T : n T (x)f'(u(x,t)) < 0} . ^ 

r 0 = r \ r, , v 

i. e., T/ corresponds to ingoing characteristics and To to outgoing ones. 

Proposition 5.1 Suppose that there is a smooth solution u(x,t ) of the initial-boundary 
value problem (37) for t > 0 and that the flux functions fi satisfy condition (39). Then 
u(x,t) satisfies the maximum norm estimate 

|W-,()IU<IMU + IWIr,.„, 

where = sup, ierj |^(x,i)|. 


Proof: 

As in the one-dimensional case, the canonical splitting can be applied repeatedly to each 
separate flux vector (/,) X| to yield 
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where 

/f ’- 11 = f 

Jo 

Due to condition (39) it follows that 

sgn (v T f'(u)) = sgn (jf v T f {v)v 2p ~ l dv'j 

for all unit vectors KR" and p = 1, Consequently, the cone condition holds for each 

p. Integration over 0 implies 

^ll«ll? = ~ 2 pJ t J 0 v T f'(v)v 2p ~ l dvds 
by the divergence theorem. Hence, 

J t W u \\lp - ~ 2 p J r J 0 v T f\v)v 2p - l dvds 

since the cone condition holds for each p. Along V j boundary data ?/>( x , f) is given, i. e., 

A r r\^\ 

j t \\ u \\ 2 p < 2 P J r j o \v T f'(v)\v 2p -'dvds. 

Let t be fixed. Define 

M = max \v T f(v )\ . 

0<V<\lp(x t T)\ 

0<T 

x6 T; 

Integration with respect to time results in 


IK, Ollg < IMS + M f' j Ti M 2r dsdr . 

Thus 

where (T/ 1 denotes the surface measure of T/. Taking the 2pth root of this inequality 
yields 

MsOllap ^ IMl2 P + (|r/|A/<)^ llV’lk.oo- 

Letting p — ♦ oo proves the proposition. □ 

Next, we consider (37) where u,/,(u) € R**. It will be assumed that the flux vectors can 
be symmetrized by a change of variables u = u(v), i. e., u'(v) is symmetric positive definite, 
and f-(v), i — are symmetric (fi(v) = /,(u(u))). The boundary conditions are 

defined in terms of the new variable v by eqs. (23), (24) (with u replaced by v = n(x,t)) 
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at every boundary boundary point x E T. The inflow and outflow boundaries T; and To 
are defined as 


T/ = {x € r : Aj(u(v(x, t))) < 0 for some j, 1 < j < d) / 4 g\ 

r G = r \ r, 

where the A ; are the eigenvalues of 

f{v,v) = ( 44 ) 

1 = 1 

i. e., A = Q T f'Q\ Q is composed of the (orthogonal) eigenvectors of /'. 

Proposition 5.2 Suppose the initial-boundary value problem (37) can be symmetrized by 
a change of variables u = u(u) such that the flux vectors /(u(u)),- = fi(v) satisfy f-v = pf, 
for some p > 0. Then any smooth solution u(x,t) satisfies a generalized energy estimate 

IW-.OIft < \M\l + f m-,r)\\l,,,dr , 

where 

ll^llr,.A/ = L 'I’Jkii’i- 

JV i 


Proof: 

Let U(v), Fi(v) E R d be solutions of 

U'v = -U + u F-v = —Fi + fi 
u(0) = 0 Fi{ 0) = f,( 0) . 

By means of the canonical splitting we get 

(v T Cr + = o. 


Since fi satisfies f-v = pf, it follows that 


F' = /' 

* p+l J% 


and thus 


(v 


T U'v)t + — — - ^2{v T fi v )x, — 0 • 


The divergence theorem yields 

d_ 

dt 


v\\f, — — f v T f'(v,v)vds = / (u>JA/u>j — UqAou>o) ds . 

u p + 1 Jr p + l Jr v 7 
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According to the definition of T / and To it follows that A/ = 0, A 0 = A > 0 on r 0 . Thus 

^IMI (j — ~ + | J r (uj — UqAouo) ds . 

The boundary conditions on T/ imply 

Integration with respect to time yields the desired estimate. □ 


6 Weak Solutions 


In all of the previous sections we have assumed the solutions to be smooth. This as- 
sumption allowed us to introduce the canonical splitting, which is a crucial step in the 
derivation of the energy estimate. Suppose now that u is a weak solution of 

«< + £ (/.)*, = o 
1 = 1 

V = Pv + Qj4> 
u(x,0) = y?(x), 


x 6 ft C R n t > 0 

x e r ; 


(45) 


where v = v(u) € IT* is a change of variables that symmetrizes eq. (45); T/ is defined 
by eq. (43), and the projection P is given by eq. (24). In general, one cannot hope for 
an energy estimate to hold since the canonical splitting no longer exists. We therefore 
assume that eq. (45) can be regularized in the following sense. Let u M = u(u M ) be a smooth 
solution of 

n n 

Ut + (/*)x, = ^ u *<*. x € ft C R n t > 0 

/j.u u + (/ — P)v = Qixjj xGT ^6) 

u(x,0) = y>(x), 

where fx > 0; u u denotes the normal derivative of u at T. 


Definition 6.1 A weak solution of the symmetrizable system (45) is called a viscous limit 

If 

lim u M (x, t) = u(x,/) (dx)-a. e., x £ 0, t > 0 

lim u M (x,f) = u(x,<) (ds x dt)-a. e., (x,t) £ T x (0, oo), 

and if these limits are bounded on ft and T x (0,f) for each fixed t; dx is the n-dimensional 
Lebesgue measure, and ds is the (n — \) -dimensional surface measure on T induced by dx; 
dt is the 1 -dimensional Lebesgue measure. 
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Proposition 6.1 Let u be a weak solution of the symmetrizable system of conservation 
laws (45)- If u is the viscous limit of solutions of (46), then u satisfies the entropy 
condition 

U(u) t + j2F,(u) x ,<0 

1 = 1 

in the integral sense; U{u) and F,{u) are given by 

U(u) = / v T u'{6v)6vd6 Ffu) = f v T f[{9v)9vd6 . 

Jo Jo 


Proof: 

We mimic the proof of Theorem 5.6 in [6]. The regularized solution satisfies 






1 = 1 


1 = 1 


where = u(n M ), v M ff = /,(u M ) = /,(u M ). Since is a smooth solution we 

can split the hyperbolic part canonically, which after scalar multiplication by leads to 


1=1 


1=1 


where 


U(v„)= f u{9v^)d9, Ffv^) = f f i (9v li )d9 . 

Jo Jo 

Define the scalar functions 

U(u M ) = vlu'iv^v^ , P.(u M ) = • 

Hence, (/'(«") = »• e -> where is 

a convex function. We thus get 


CK), + v = r E C(u%.r. -*£,<? C'VK • 

t=i «=l *=i 

Since U is convex it follows that u£ r (/"(u M )u£ i > 0, that is, 

£/(»')< + E < /* E c(u B )„,, . 

1=1 1=1 

Integration over f l x (0,t) and the divergence theorem imply 

J U(u^)dx-J U{<p)dx + J j^v T Fiu^dsdr < p J^U^u^dsdr , (47) 
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where U„ denotes the outward normal derivative, and 


»=i 


But 


»jj J r U y {u»)dsdT = nf Q J r U'^Kdsdr . 

The normal derivative u £ can be eliminated by means of the boundary conditions. Thus 

fi £ j r Uy(u M )dsdr = £ J vliQfy -{I- P^v^dsdr . 

Using definition (6.1) we can apply the Lebesgue dominated convergence theorem as 

H —* 0 

lim n J o U v {u?)dsdT = jf J v T (Q I tl> - (I - P)v)dsdr . 

At the outflow boundary To we have Qj = 0, P = /, i. e., 

^ Jo Jr Uv{ u ^)dsd T = J jf — [I — P)v)dsdr . 

On T/, however, the boundary conditions v = Pu + Qjty are satisfied ( dt x ds)- a. e. 
Consequently, 

lim/x / / = 0 . 

m—o Jo Jr 

The left hand side of eq. (47) can be handled analogously. We have thus arrived at 

j Q U{u)dx - jf U(ip)dx + fj J v T F(u)dsdr < 0 , 

which concludes the proof. □ 

As a consequence of proposition (6.1) we have 

Proposition 6.2 Let u be a weak solution of the symmetrizable system of conservation 
laws (45) where the flux vectors satisfy f'v = pf for some p > 0, i = 1 , . . . , n. If u is the 
viscous limit of solutions of (46), then v = v(u) satisfies a generalized energy estimate 

IM'.OI6<IM&+ ? | T X‘lW-,r)|g 1-Al rfr. 


Proof: 

From proposition (6.1) it follows that 

U(u)dx — J U(<p)dx < — u v T F{u)dsdr , 
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where 


U(u) = r v T u(0v)0vd0 = v T U'v , U(v) = f u(0v)d0 

Jo Jo 

v T F(u) = ViFi{u ) = u i I v T f' i {0v)0vd0 . 

.=1 «=i Jo 

Using f-v = pfi, the identity 

f « > T fl(0v)0vd$ = -L-v T fl(v)v 

Jo p + 1 

follows easily by integration by parts. Thus 

u t F(u) = — J— ^ = — J-7 v T f'{v, u)v . 

p + I P + * 

The entropy condition can then be expressed as 

IM-. Oiler < M \ 2 V~y+\l J r V T f'(v,v)vdsdT. 

The remainder of the proof is identical to that of proposition (5.2). O 

The essential step in the proof of proposition (6.1) was the establishment of scalar 
functions U{u), F,(u) ( U(u ) convex) satisfying U tT f i = F? . In the scalar case we know 
that there is a sequence {U p , F^}, p = 1,... such that U' p f' t = F t ^\ U p and F/ p ' are 
given by 

U,(u) = J-« 2r , Fj\ u) = f f'lvy-'dv . (48) 

2 p Jo 

For each pair {F p ,F/ p) } proposition (6.1) is true, i. e., 

UM t + ±F?\uU<0 (49) 

t = l 

holds for all weak solutions u of eq. (37) that are viscous limits of the regularized problem. 
If condition (39) is met, it follows that the cone condition (5.1) is valid for all p. Equations 
(48), (49) thus imply an energy estimate for each p. Taking the limit as p — > oc yields 

Proposition 6.3 Let u € R be a weak solution of the scalar conservation law (37), where 
the flux functions /, satisfy condition (39). If u is the viscous limit of the regularized 
problem, then u(x,t) satisfies the maximum norm estimate 

IM-.OIloO < Moo + IMIr/.OO • 
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Thus far, it has been tacitly assumed that the sign of the eigenvalues Aj(u(x,<)) is 
independent of t for each i G T. This implies that for a given boundary point the rank 
Qi ( Qo ) is constant in time, which means that the rank of P also is constant in time 
(recall that the coupling operator S has to be w small”). As a consequence, the inflow 
boundary T/ is time independent. 

The estimates obtained in the previous sections have been derived from an extra con- 
servation law (or possibly a sequence of conservation laws) subject to the same boundary 
conditions as the original one. To study the impact of piecewise rank constant boundary 
conditions it suffices to consider eq. (45). Assume that the rank of P is constant on [0, fj) 
and [<i,oo). According to proposition (6.2) we have 


\W-,t)\\b<M\o + -^ J o HV> ( 0 ) (-,r)||^ 0 ) ^ 0 ) dr, 0 <<<<!, 


and 


ll u (') Oil# ^ IM'j i £ llV ,(1) (-,r)||^ J) A (,)dr, ti<t< 


oo 


Now, since v — v(u) is the viscous limit of a regularized problem it follows immediately 
that 


which yields 


Define 


lim|H-,<)ll& = IH-,*i)lli/, 

f v (* T ^i) I lt> < IMIgr + ~ ~ j j o l|0 (O) (-,7')||p(o) A (o)dr . 

r _/ri 0) o <t<t x 

;_ l «,<* 


with similar definitions of ip and A/. Thus 


IN-'OIIJ < Mo + “i jf M;r)\\ 2 r lAl dr, t >0, 

which is identical to the estimate of proposition (6.2). The maximum norm estimates 
follow from a sequence of inequalities of the above type. Consequently, all energy estimates 
and maximum norm estimates hold when the boundary conditions have piecewise constant 
rank in time. 


7 Discussion 


Our aim has been to derive energy estimates for the initial-boundary value problem for 
nonlinear conservation laws by means of a few basic principles: the canonical splitting 
and the cone condition. These principles apply to scalar conservation laws and systems, 
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to one or more space dimensions (even the time dimension), to smooth solutions and 
weak solutions that are viscous limits of vanishing viscosity solutions. For symmetrizable 
systems the canonical splitting gave rise to a new scalar conservation law 

U(u) t + j2F t (u) Xi = (<) 0 (50) 

t=i 

where the convex function U(u) and F t (u) are related via U fT f- = F fi{ u ) 1S the flux of 
the original conservation law. Thus, the canonical splitting implies an entropy condition. 
Furthermore, we have obtained an explicit expression for the entropy function U(u) and 
the entropy flux F t (u). 

For the initial value problem a generalized energy estimate follows immediately by 
integration of eq. (50), since 

U(u) = v T U f (v)v , 

where U f (v) is symmetric positive definite; v = u(u) is the change of variables that 
symmetrizes the original conservation law. If we want an energy estimate for the initial 
boundary problem, then we need to impose a cone condition on the flux vectors /*(«). 
The cone condition ensures that potentially dangerous boundary terms will be eliminated 
by enforcing the boundary conditions. The cone conditions will, in general, impose rather 
stringent conditions on /j(u), even in the scalar case. Convexity, however, is not needed. 
This contrasts with the initial value problem, in which case one can prove existence and 
maximum norm estimates for single conservation laws under very general hypotheses (/(«) 
differentiable, initial data ip(x) € I^R") H L~(R n ) H £K(R n ))[2]. Albeit a restrictive con- 
dition, the cone condition is fulfilled for many scalar equations that occur in practice, 
such as the Burger’s equation and the Buckley-Leverett equation. For systems a partic- 
ularly important class of flux vectors is given by those that are positively homogeneous 
functions of u, i. e., f'v = pf for some p > 0. The Euler equations belong to this category. 
Summing up, if a cone condition holds, then a generalized energy estimate follows from 
eq. (50) and the characteristic boundary conditions. Thus, the energy estimate appears 
as a special case of the entropy condition (50). In particular, there exists a generalized 
energy estimate iff the system is symmetrizable and a cone condition holds. 

We saw for a single conservation law that it was possible to generate a sequence 
of entropy conditions from which we could deduce a maximum norm estimate. It is 
an intriguing question whether systems of conservation laws that originate from physics 
possess a sequence of entropy conditions such that a maximum norm estimate would 
follow. As mentioned earlier, Lax [5] has constructed such a sequence for 2 x 2-systems. 
Harten [3] has introduced a family of entropies for the Euler equations such that symmetry 
and homogeneity are preserved. It remains to be shown if a family can be found that 
satisfies the necessary growth conditions. Finally, we conclude this discussion by noting 
that many of the energy estimates carry over to the semi-discrete case by extending the 
technique in [9] to the nonlinear case. This is the topic of a forthcoming paper. 
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